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OF LONG-RANGE BALLISTIC MISSILES 
By A. J. Eggers, Jr- 

SUMMARY 


It is demonstrated with the aid of si mi litude arguments that a model 
launched from a hypervelocity gun upstream through a special supersonic 
nozzle should experience aerodynamic heating and resulting thermal stresses 
like those encountered by a long-range ballistic missile entering the 
earth’s atmosphere. This demonstration hinges on the requirements that 
model and missile be geometrically similar and made of the same material, 
and that they have the same flight speed and Reynolds number (based on 
conditions jiist outside the boundary layer) at corresponding points in 
their trajectories. The hypeirvelocity gun provides the model with the 
required initial speed, while the nozzle scales the atmosphere, in terms 
of density variation, to provide the model with required speeds and 
Reynolds numbers over its entire trajectory. Since both the motion and 
aerodynamic heating of a missile tend to be simulated in the model tests, 
this combination of hypervelocity gun and supersonic nozzle is termed an 
atmospheric entry simulator. 


INTRODUCTION 


The aerodynamic heating of a long-range ballistic missile entering 
the earth’s atmosphere poses problems of so serious a magnitude that the 
success or failure of the missile may well hinge on their solution. The 
extent to which these problems are solved can be fully determined only by 
flight tests, but such tests are, per se, very difficult, time consuming, 
and costly. It is appropriate to inquire, therefore, if a method can be 
devised for simulating with relatively simple equipment on the groimd, 
the aerodynamic heating and resulting thermal stresses in a bsillistic ' ^ 
missile. This question forms the starting point of the present paper 
which undertakes first to establish conditions of simulation and then to 
determine a practical method of simulation. In the course of events it 
will be discovered that simulation of heating goes hand in hand with simu- 
lation of motion of a missile; accordingly, the subject of this paper is, 
more generally, an atmospheric entry simulator. 
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SYMBOLS 

reference area for drag evaliiatlon, sq. ft 
rate constant (see eq. (Al)), sec"^ 
rate constant (see eq. (A2)), cu ft/sec 
drag coefficient, dimensionless . 

local skin-friction coefficient based on conditions just 
outside of boimdary layer, dimensionless 


equivalent skin-friction coefficient, 
dimensionless 



typical model or missile dimension, ft 

Young's modTolus, Ib/sq ft 

heat transferred per unit area, ft-lb/ft^ 

diffusivity, sq ft/sec — ' 

thermal conductivity, ft-lb/sec ft^(°R/ft) 

length of nozzle, ft 

flight tfech number, dimensionless 

mass of missile or model, slugs ' 

concentration of particles (i.e., number of particles per 
unit volume), ft”® 

distance normal to boxmdary, ft 

reduced normal distance, n/D, dimensionless 

nozzle reservoir pressiu’e, Ib/sq ft 

total heat transferred, ft-lb 

radius of ciirvature of body svirface at stagnation point, ft 


r 
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surface area, sq ft 
temperature , °R 

nozzle reservoir temperatirre , °F 
time , sec 

reduced time, sec/sq ft 
velocity, ft/sec 
altitude, ft 

rectangular coordinates, ft (except when appearing as sub- 
scripts on stresses) 

reduced coordinates (x = x/d, y = y/D, z = z/d), dimensionless 

coefficient of thermal expansion, ft/ft-9R 

constant in density-altitude relation, ft”^(see eq. (l) ) 

typical flow deflection angle, deg 

angle of flight path with respect to horizontal, deg 

Poisson's ratio, dimensionless 

air density, slugs /cu ft 

reference density, slugs/cu ft (reservoir density in nozzle) 
tensile or ccanpresslve stresses, Ib/sq ft (see sketch, p. 9) 
shear stresses, Ib/sq ft (see sketch, p. 9, and note 

Txy = Tyx^ Txz = '''zX^ '''yz ^ '’’zy) 

direction cosines of normal to a boundary, dimensionless 


Subscripts 


av average conditions 

e conditions at "entrance" to earth’s atmosphere or simulator 


I 


local conditions just outside of boundary layer 



k 
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mi missile 

mo model 

s stagnation conditions 


ANALYSIS 


In order to clsLrify the basic requirements of a simulator, this study 
is initiated with a review of some of the more important causes and effects 
of aerodynamic heating of long-range ballistic missiles. Following this 
review, the conditions of simulation will be set forth in a mathematicaJ. 
form designed to facilitate the choice of a practical simulator. 


Motion and Heating of Long-Range Ballistic Missiles 


It is a fundamental characteristic of ballistic missiles (see ref. l) 
that speed begets range, with the result that hypervelocities in excess 
of 10,000 feet per second are required in order to obtain long ranges in 
excess of 1,000 miles. A long-range ballistic missile first attains 
iiyper velocities near the end of powered flight. In this phase of flight 
and throughout the large majority of unpowered flight the vehicle shoiild 
pass more or less imimpeded through the rarefied upper atmosphere of the 
earth, corresponding to altitudes in excess of several hundred thousand 
feet.^ Its trajectory terminates, however, with a very rapid descent 
throxigh all or part of the earth's relatively dense lower atmosphere. In 
this phase of flight, termed the atmospheric entry, retardation and severe 
aerodynamic heating of the missile can almost- certainly be expected to 
occur (see ref. 2). 

Retardation during atmospheric entry is caused by the combined action 
of pressure and viscous forces, while aerodynamic heating stems in the 
main from work done by viscous forces. In both cases it is aerodynamic 
rather than gravity forces which play the predominant role, with the result 
that motion and heating of the missile emerge as closely related phenomena. 
Thus changes in missile shape which affect motion will also affect heating. 
This fact can, as discussed, in reference 2, be exploited in the design of 
missiles with reduced aerodynamic heating. 1116 potential for excessive 
heating remains, however, as an -unavoidable property of the long-range bal- 
listic missile which enters the earth’s atmosphere at hypervelocity. 


^It is presumed that the missile is, unlike the usual meteor, so 
large (say of the order of feet in typical dimension) that free molecule 
phenomena play a minor role in its motion and heating. 


NA.CA RM A55H5 


5 


Excessive heating can have several effects. First and perhaps most 
serious of these effects is the development of high thermal stresses in 
the structure of the missile. Hhese stresses tend, for example, to far 
overshadow the pressure induced stresses. In addition, there is the 
natural weakening of the missile material at high temperatiares , so that 
structural failure may occur during atmospheric entry. There is, of 
course, the further possibility that intense aerodynamic heating will, 
as with meteors, cause burning and the ultimate destruction of the missile. 

These, then, are some of the important causes and effects of aero- 
dynamic heating of a long-range ballistic missile entering the atmosphere. 
They suggest that an atmospheric entry simulator might logically include 
means of simulating missile velocity, missile configuration (e.g., shape), 
and the lower portion of the earth's atmosphere. We will proceed from 
this suggestion to formulate om* conditions of simulation. 


Conditions of Simulation 


The analysis to follow tacitly pres\mies the validity of many time- 
tested assimptions of aerodynamics, themodynamlcs, and solid mechanics. 

In addition, however, it is predicated on the assumptions that during 
atmospheric entry (l) radiation has a secondary effect on missile heating, 
( 2 ) gravity has a secondary effect on missile motion, ( 3 ) the flow-field 
freeze principle of Oswatltsch (ref. 3) holds for the i^ssile, and the 
thermal properties (e.g., specific heats) and transport properties (i.e., 
thermal conductivity and viscosity) of air are fvinctlons of temperatiire 
only. The first two assumptions are suggested by the calculations of 
references 1 and 2,^ Tihile the third assumption hinges essentially on the 
requirement that the square of the hypersonic similarity parameter for a 
missile be large compared to 1 d-uring entry (i.e., M^sln^b > > l) . For 
missiles of normal slenderness, the hyperspeed of entry tends to Insure 
the satisfying of this requirement with the result that flight Mach num- 
ber loses its significance as an Important similarity parameter.® The 
last assumption is restrictive only in the event air tends to dissociate 


^The first assiunption is, of coiurse, best suited to "relatively light" 
missiles which are designed on the so-called "heat sink" principle, or, 
more generally, which are designed to maintain relatively cool surfaces. 

®It is demonstrated in reference 3 that Miach number and stream angle 
in the disturbed flow are independent of flight Mach number, provided 
M®sin®8 > > 1, and provided the air behaves ideally. In the event non- 
ideal behavior, like charges in the specific heats, occurs as a result of 
high temperat'ures in the dlstiirbed flow, then flight velocity tends to 
replace flight Mach number as the important index of motion since it is 
thrcu^ this velocity (via kinetic energy) that hi^ disturbed air tem- 
pera tiures are created. 
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(or possibly ionize) and it will therefore be treated in a disciission 
of these phenomena later in the paper (see section on "Performance 
Limitations" ) . 

Now it is convenient in discussing similitude to imagine a model 
coinaterpart to the missile and a test chamber counterp^t to the atmos- 
phere. Furthermore it is permissible for our purposes to proceed from 
the simplified equations of reference 2 for the convective heating of 
ballistic missiles. Thus for an isothermal atmosphere (which closely 
approximates the earth’s lower atmosphere, see ref. 2) 


Po 



and there follows : 

(a) Heat absorbed per unit mass at altitude y* 

- = f (%r) 

where 


= 


Ve^e 


CpPpA 
|3m sin 0g 


,-Py’ 


( 2 ) 


(3) 


(l) Average rate of heat transfer per unit area 


dt 


Cf 'p 


o^e' 


4 



3 CpPpA _pyl 

^ ^m sin &Q 


(^) 


(c) 


Rate of heat transfer to stagnation point 


dH I - ^ _ 3 CpPpA py. 

= 6.8xl0"®y^ Ve®e ^ e ^ (5) 

According to equation (l) the density of the air in the test chamber 
must vary exponentially with the distance corresponding to altitude in the 
atmosphere . The more general implication is , of course, that the test 
chamber must duplicate variations of p/p^ in the alanosphere, whatever 
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they may be, although the absolute magnitudes of p and Pq may be quite 
different from those in the atmosphere- The static temperature of the 
air in the test chamber is, as in the case of the atmosphere, presumed to 
be small by comparison to missile recovery temperature (see ref. 2 in 
connection with this point as it relates to the derivation of eqs. (2) 
through ( 5 ) ) • 

It will be stipulated now that model and missile be geometrically 
similar in structure and configuration, and made of the same material. 
Furthermore, the condition is imposed that the model enter the test cham- 
ber at the same speed and temperature as the missile enters the aianosphere . 
Finally, it is required that model and missile have the same Reynolds num- 
bers (based on local conditions outside the boundary layer) at correspond- 
ing points Py' in their trajectories. By corresponding points it is 
meant where the product Py' is the same for model and missile. It should 
be recognized, of course, that, in general, p and y' wi3J. individually be 
grossly different for model and missile. 

It follows frcm these requirements and eqviation (2) that the heat 
transfer per unit mass Q/m to model and missile will be the same at 
corresponding points Py' provided V is the same, since Cf ' and 
Cf'S/CoA are the same. But from equation ( 3 ) the velocity V will be 
the same at corresponding Py* only if Cj)PQA/pm sin 0^ is the same. 

If the subscript mo refers to model and mi to missile, then the last 
provision may be written 


(PoD)mo (y’/sin 9^)^ 


But model and missile Reynolds numbers, velocities, and disturbed air tem- 
peratures^ are the same, hence 


This observation with regard to disturbed air temperatures is easily 
verified by considering flow near the surface of a missile with a stagna- 
tion point. Assuming for simplicity that air in the disturbed flow 
behaves ideally, we have at the stagnation point Tg k V^/2Cp since 

> > 1. Then the temperat\ire of the air just outside the boundary layer 

is given by the expression Tj = — — (l + ^Ij • Hence if V and Mj 

2Cp \ 2 / 


are the same for model and missile, then Tj and Vj are the 
pendent of ambient air temperature. We are assured of equal 
freeze principle (see footnote 


same, 

Mj's 


inde- 
by the 
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and equation (6) may be -written 



Elis expression fixes the length L of -bhe test chamber in terms of the 
portion of the a-tmosphere to be simulated therein and the ratio of model 
to missile size.^ If equation (8) and the previously set forth require- 
ments are satisfied, then model and missile should experience equal heat 
transfer per unit mass, and hence equal average ten^jerature rise at corre- 
sponding points in their trajectories. These quantities are significant, 
of ■course, because they tend to determine -whether a missile -will melt or 
perhaps burn during flight. 

The next question is how do tile heat-transfer rates compare in the 
case of model and missile? It is easily deduced from equations (i)-) a.-nfl ( 5 ) 
and the conditions for eq-ual heat transfer per unit mass that 


and 



Aal 

%o 





Dmo Vdt 


( 10 ) 


at corresponding points Py'. That is to say, the a-verage and stagnation 
point heat-transfer rates are higher for the model in proportion to the 
ratio of missile to model size . But perhaps the foremost Importance of 
heat-transfer rates is, as discussed earlier, in how they Influence 
thenial stresses in the missile structure and, for example, lead to abla- 
tion of surface material. Evidently, then, it would be most desirable 
if equations ( 9 ) and (lO) implied equal theimal stresses in model and 


®Since model and missile velocities are the same at corresponding 
points in their trajectories, it follows from equation (8) that the time 
of flight in the test chamber is reduced belo-W that of a-tmospheric entry 
by the ratio (Dmo/Anl)^* “ " 
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missile. Eiis possibility is easily checked using modified equilibrivim 
thermal-stress equations for an unrestrained isotropic elastic body (see 
ref . (4) and sketch) 


2 



Stress notation 


Thus 



Sr^y 


^Txz 

(xE BT Q 

Bx 


T 

,dz 

1 - 2v Sx 

aCTy 

9rxy 

J. 

^Ty2 

oE St ^ 

Sy 

8x 


5z 

1 - 2v Sy 

- -1 

^Txz 

_L 

STyz 

oE Q 

5z 

3x 

"T 

Sy 

1 - 2v 5z J 


( 11 ) 


where the stresses are not true thermal stresses, but rather they are the 

stresses produced by the "body forces," etc. 

1 - 2v Sx 
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!nhe corresponding compatibility conditions are 


where 


2 iSlr oE/v o 

^ 1 + v l- 2 v\l-v 

+ 2 aff 

1 - 2V V- “ b^y 




+ 


V^Tyz + 


V^-TXZ + 


+ 


1 b^lf 

oE / V 

1 + v 

1-2V \1- V .. 

1 ^ 

^ oE 

l+v SySz 

1 - 2v Sy^z 

1 _ 

oE b^T 

l+v Sx5z 

1 - 2v Sxdz 

1 

^ oE 

l+v 3x3y 

1 - 2v bxbf 




bx^ by^ bz^ 


i 


Sz 


-2 




> (12a) 


and the boimdary conditions are 


■*■ ’’’xv^ ■*■ ''xz^ “ 


ctET 

xyt ■ -xzi* “ 1 _ 2v 


oET 


■*■ "^yz^ ■*■ "^xyS ^ Y^2v 
a,t * + xy,, = ^ 5 




(12b) 
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The heat-flow equation is (see ref. ( 5 )) 


by^ K ^t 


( 13 ) 


with boundary conditions 


and (probably) 


(^wall) 


initial 


known 



outer 

surface 




= O 

inne r 

surface 




> 

> 


(ll^) 


Now the thermal stresses are given by superposing the "hydrostatic" pres- 
sure otET/l-2v on solutions to equations (ll) consistent with the com- 
patibility and boundary conditions (eq_s. (12)) and the solution to the 
heat-flow equation ( 13 ) along with its boundary conditions (eqs. l4) . 

But if equations (9) and (lO) along with the requirements necessary for 
their develojanent are satisfied, then equations (ll) through (l^) (and 
the "hydrostatic" pressure) are mathematically identical for model and 
missile. Hence, the thermal stresses must be the same in model and 
missile at corresponding Py’ in their trajectories. 

Necessary conditions for simulating aerodynamic heating and resulting 
thermal stresses in a ballistic missile are then, according to this anal- 
ysis, as follows. First, model and missile must be geometrically similar 
and made of the same material. In addition, they must have the same 
flight speeds and Reynolds nximbers (based on local conditions outside the 
boundary layer) at corresponding points in their trajectories. Finally, 
in order to meet these conditions and to insirre equal heating of model 
and missile, the test chamber must contain air at relatively low tempera- 
ture and with variations in p/Pq along the model fli^t path equal to 
those in the atmosphere along the missile flight path. 

With this knowledge we are in a position to consider the practical 
problems of simulating atmospheric entry. 
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A PRACTICAL ATMOSPHERIC ENTRY SIMULATOR 


It is appropriate to determine first how to- provide a model with 
the correct Initial hypervelocity required for simulation. For this pur- 
pose It Is suggested that a hypervelocity gun can be employed. A. C. 
Charters and his group at Ames Laboratory (see ref. 6) have, for example, 
launched a 22-callber model from a hypervelocity (llght-gas) gun at speeds 
up to 15,^00 feet per second. This speed Is 60 percent of satellite 
speed and corresponds to that of a ballistic missile (see ref. l) with up 
to 2,000 miles range. Substantially higher speeds (corresponding to 
longer ranges) appear to be realizable with light-gas guns. Both speed 
and lange fall, then, into the categories of Interest in this paper. 

The next question is how to provide a model test chamber which simu- 
lates the lower portion of the earth’s atmosphere. For this purpose it 
is suggested that a special supersonic nozzle can be used to advantage. 

To illustrate, it was foimd in reference 2 (see fig. 4 therein) that the 
major part of the aerodynamic heating of a ballistic missile entering 
the atmosphere occurs over a range, of about 100,000 feet in altitude. 

In this altitude range p/pq varies by a factor of about 10“^. A corre- 
sponding variation in density can be obtained between the settling cham- 
ber and exit section of a Mach number 5 supersonic nozzle (see ref. 7) • 
Imagine, then, a hypervelocity gun positioned to launch a model upstream 
along the axis of such a nozzle. The nozzle is designed to provide an 
essentially exponential variation in density along its axis, the density 
decreasing from a very high value in the settling chamber to a very low 
vaJ-ue at the exit. Accordingly, a model proceeding upstream through the 
nozzle encounters an increasingly dense atmosphere like that presented 
by the earth to a descending ballistic missile. Now, to be sure, unlike 
atmospheric air, the air in the nozzle is in motion. However, the air 
velocity is small by comparison to the hypervelocity of the model, and 
therefore this difference between nozzle air and atmospheric air should 
mar only slightly the function of the nozzle as a test chamber.® 

The combination’ of hypervelocity gun and supersonic nozzle merits 
attention, then, as an atmospheric entry simulator. An example simulator 
of this type will therefore be considered next. 


®It is interesting to note that the supersonic nozzle when used in 
the proposed manner has compensating features ; namely, as model velocity 
decreases , the air velocity decreases and, as model recovery temperature 
decreases, the air static temperature rises up toward atmospheric air 
temperature. These features are, of course, favorable to simulating 
atmospheric entry in that they tend to preserve the required similarity 
between atmospheric air and nozzle air. In this manner, too, the pos- 
sibility of correcting for small differences between these media is 
enhanced . 


(lODtBLDBNTlAE^ 
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Example Simulator 


Consistent with the previous discussion^ a Mach number 5 nozzle is 
chosen to simulate the earth's atmosphere over 100,000 feet of altitude. 
Now- let us assume that ballistic missiles with up to ^,000 miles range 
are to be studied with this simulator. The corresponding rajige of atmos- 
pheric entran.ce angles 6 q is from down to about 30° (see ref. l) . 

The length of the nozzle is fixed according to equation (8) by the Trw.viTmuTi 
values of (yg'/sin ^eVi l|moAki* value of (yg'/sin 9e)ijii 

from the above specifications, ^ = 200,000 feet. The maximum value 

sin 30° 

of Dmn/Pmi will be dictated by the size of the largest model which can 
be launched by the special gun available and by the size of the smallest 
missile to be simulated. For the purpose of this discussion it suffices 
to observe that the model size will probably be of the order of a fraction 
of an inch, while the missile size will probably be of the order of sev- 
eral feet. It follows that the maximum value of Dmo /Dmi should be of 
the order of 10“^. In this event we have from equation (8) that the 
length of the nozzle is of the order of 20 feet, and from equation (7) 
the nozzle stagnation density is of the order of 100 times sea level air 
density. A missile 3 feet in diameter and weighing 5^000 pounds would 
thus be simulated by a model O .38 inch in diameter and weighing 0.005 
pound. 

Cta. the basis of these considerations our example simulator might 
appear seme thing like the one shown in figure 1. The required nozzle 
stagnation densities ane obtained for running times of the order of a 
second by using a settling chamber charged with high pressure air. A 
wide range of settling-chamber pressures and, hence, densities makes it 
possible to vaiy the range of altitudes simulated. The storage heater 
in the settling chamber maintains the air temperatures above those for 
which liquefaction can occur in the nozzle. The nozzle contour will tend 
to have very small slopes in the streamwise direction. Hence it is 
anticipated that flow in the nozzle wi3J. not depart radically from the 
one-dlme 2 islonal type. Air from the nozzle passes into a vacuum tank 
which is of sufficient size to maintain the nozzle compression ratios 
required for supersonic flow d-urlng the course of a test. Observation 
windows ane located at short Intervals along the nozzle side walls to 
permit photographing models in fll^t and detemnining effects of aero- 
dynamic heating. In this regard it would quite likely prove desirable 
to employ a spectrograph (or spectrographs) to identify the sources of 
radiation energy emitted in the vicinity of models. The simulator would 
be instrumented for meastiring settlisag-chamber pressinre and temperature 
and, of course, the time-distance history of models. Models would be 
la u nched from a hypervelocity gun located at or near the end of the noz- 
zle, and they would be arrested in a catcher located in the settling 
chamber . 
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We have, then, scone idea of how the proposed simulator might be 
employed in practice. It is important, however, to be aware of the 
points of difficulty which may limit the performance of the device as 
a simulator. 


Performance Limitations 


First, it shouild be noted that pressure -induced stresses in the model 
will be higher than those in the missile by the ratio Dmi /^mn • ^Chis 
point may prove troublesome, although not unduly so if these stresses are 
very low (as they tend to be) in the missile.' It should be remembered 
too that by proper model design (e.g., the use of internal pressurization 
or other bracing) this problem can be minimized. Careful attention must, 
of course, be given to both model and sabot design from the standpoint of 
minimizing stresses in the model during launching. 

Another point of difficulty may be encountered if material properties 
(e.g., yield point) are a significant function of time under rapid heating 
conditions. There is some indication (see ref. 8 and papers cited therein) 
that at the very high heating rates of long-range ballistic missiles 
(corresponding to temperatxire-rise rates of the order of 100° F/sec and 
more), the Importance of time is small. It is indicated too that this 
remains true, and, more important, that material properties remain essen- 
tially the same at the much higher heating rates produced in the simulator.'^ 
Far more information is needed, however, befora the significance of time 
in the sense of this discussion can be fially assessed. 

A further question which should be raised is whether or not a body 
can distort fast enough with changes in temperature to remain in stress - 
strain equilibrium (i.e., acceleration terms negligible) as was assumed 
in deriving equations (ll) and (l2) . This situation wo\ild be more serious 
in the case of the model than in the case of the missile. A crude check 
on the matter is easily obtained for a body wMch behaves elastically. 

To illustrate, the time required for the model to make small adjustments 
in stress-strain equilibrium should be of the same order as the time 
required for an elastic wave to travel the length of the model. For a 
typical steel model (D=1/3 in.) in the simulator, this latter time would 
be about 10“ ® seconds . If we assime that the model experiences a total 
temperature- rise of 1500° F while traveling through the simulator, then 
10“® seconds is also the time required for only about 1° temperature 
change in the model material. Accordingly, only small adjustments in 
stress-strain equilibrium are evidently required in this time, and hence 


Y 

This discussion shoiild not be construedT to mean that material prop- 
erties are the same under conditions of high rates of heating as they 
are under steady state conditions (see ref. 8). 
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equilibrium should tend to be realized in model as well as missile. If 
time -dependent plastic deformation should become significant, then this 
statement obviously no longer holds.® Thus, for exajuple, the simulator 
may not (in view of its foreshortened time scale) duplicate more than 
qualitatively a fracture process (see ref . 9) although it should tend to 
duplicate thermal deformation up to and including the beginning of 
f ractuure . 

There is, too, the possibility of a missile being aerodynamically 
heated to temperatures where it will burn during descent through the 
earth's atmosphere. The simulator should tend to duplicate conditions 
leading up to this phenomenon j however, there is a question as to how 
well buurnlng would be duplicated. The ccxnplexity and lack of complete 
understanding of metal burning (see, e.g., ref. lO) preclude the pos- 
sibility of obtaining a quantitative answer to this question at the 
present time . From the qualitative viewpoint it is reasonable to expect 
that the effects of increased partial pressure of oxygen (acting to 
Increase burning rate) and the reduced model size (acting to decrease 
the amount of material to be burned) should combine in the si mu lator to 
compensate for its foreshortened time scale, thereby more nearly provid- 
ing simulation of missile burning. Burning should be accurately simu- 
lated for cases where burning rate is directly proportional to the par- 
tial pressure of oxygen. 

As a final point, it is appropriate to consider the matter of dis- 
sociation and association of the oxygen and nitrogen in air. The simula- 
tor produces essentially the same disturbed air temperatures as the mis- 
sile entering the atmosphere. Accordingly, the potential for dissociation 
and association is duplicated by the simulator. At the present time, 
however, it is felt that these two phenomena obey different rate laws 
(see Appendix) . In this event, neither the phenomena nor their effects 
on heating a ballistic missile can be duplicated except by 1 to 1 
simulation (i.e., the equivalent of fli^t tests with the full-scale 
missile) . An indication of the possible error that this situation may 
introduce in tests with the proposed simulator can be obtained from the 
calculations of reference 11 which suggest that the net effect of equi- 
librium dissociation and association in free flight may be to increase 
the rate of heat transfer to a stagnation point by only about 50 percent.® 
This increase should be well within the design safety factor of a missile. 


®Thus with a slow process, like creep, model defoimation would not 
simulate missile defozmiation (see, e.g., ref. 9) ^ t)ut rather, it would be 
less by the ratio of (Dmn/ Pmi )^ due to the foreshortened time scale in 
the simulator. Even in the case of the missile, however, the entry time 
is so short (of the order of seconds) it seems unlikely that creep could 
play a major role in defomning the vehicle. 

^Dissociation and association may have especially marked effects 
on stagnation point heat transfer since these phenomena can strongly 
affect both the inviscld and viscous flows in this region. 
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CONCLUDING REMARKS 


It has been found that an atmospheric entry simulator consisting of 
a hypervelocity gun combined with a special supersonic nozzle may be used 
to study the conditions of extreme heat transfer and thermal stress which 
introduce such serious problems in the flight of a long-range ballistic 
missile. The effects of aerodynamic heating on the model can be observed 
with relative ease, and further, the tests can be conducted at a cost 
which is negligible by comparison to that of flight tests. Indeed, in 
the simplest test, the simulator could provide with one photograph of a 
model rather substantial evidence as to whether or not the corresponding 
missile would remain essentially intact while traversing the atmosphere. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Sept. 15, 1955 


ngp^rmMiPT.A-T-.^ 
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APPENDIX A 

DISSOCIATION AND ASSOCIATION RATES 

The dissociation of O 2 and Na is thou^t (see, e.g., refs. 12 and I 3 ) 
to obey the linear law 


^ = Ci(T)N (Al) 

where N is the concentration (or density in particles per unit volime) 
and Ci(t) is the rate constant which usually depends only upon tempera- 
ture. Cn the other hand, the corresponding association process (e.g., 

0 + 0 — > O 2 ) is thou^t to obey the second-order law^ 


^ = C2(T)N^ (A2) 

dt 

Now disturbed air velocities are the same in the case. of model and missile, 
while disturbed air densities are higher in the case of the model by the 
ratio of Dmi /Ano • It follows then that the simulator will tend to dupli- 
cate a rate process in which the percentage rate of change of concentra- 
tion of a. given type of particle is proportional to the concentration of 
that particle, namely. 


1 ^ 
N dt 


c(t)n 


(A3) 


But the simulator tends to duplicate temperatiore and type of particle; 
hence, it should duplicate the "association" rate process given by equa- 
tion (A2) since C(t) = 02 ( 1 ). It will apparently not, however, dupli- 
cate the "dissociation" rate process given by equation (Al) . 


^Actiially there is seme question as to the correctness of either 
equation (Al) or (A2) (see e.g., ref. (l4) ) for pure media and there is 
the fiirther complication of Impurities which co\ild lead, for example, 
to a third-order rate process. - 
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